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Abstract. The general exact solution of the Einstein-matter field equa¬ 
tions describing spherically symmetric shells satisfying an equation of 
state in closed form is discussed under general assumptions of physical 
reasonableness. The solutions split into two classes: a class of “astro- 
physically interesting” solutions describing “ordinary” matter with pos¬ 
itive density and pressure, and a class of “phantom-like” solutions with 
positive density but negative active gravitational mass, which can also 
be of interest in several “very strong fields” regimes. Known results on 
linear-barotropic equations of state are recovered as particular cases. 


1. Introduction 

Relativistic shells have found applications in an impressive list of differ¬ 
ent scenarios, such as sources of vacuum gravitational fields llTl-ll4l. cos¬ 
mological models Q and 0, gravitational collapse 0, canonical quan¬ 
tization of gravity |8} and @, higher dimensional models in brane world, 
quantum gravitational collapse fflOl . entropy of black holes, motion of thin 
clouds surrounding an exploding star ITU , sources of “phantom” fields and 
dark energy lH2l . In spite of the wide range of their applications, very few 
analytical solutions in closed form are known even in spherical symme¬ 
try. Exceptions are the dust (i.e. zero pressure) shells, which were com¬ 
pletely analyzed already in the pioneering work by Israel lfT3l . [fRI . and 
solutions with linear-barotropic equations of state of the form P = T 0 cr 
with a constant r 0 . In the present paper we discuss the general exact so¬ 
lution for shells satisfying an equation of state in closed form P = r(cr) 
with some general assumption of physical reasonableness on the behavior 
of the function T. Our solutions include as particular cases the linear e.o.s. 
cited above but also, for instance, non-linear polytropic equations of state 
of the form r(cr) = r 0 cr l/ lfl5l . It turns out that the qualitative behavior of 
spherical shells always presents two “branches”, an “astro-physically rele¬ 
vant” branch in which the weak energy condition is strictly satisfied, and 
a “phantom” branch in which the density remains positive but the active 
gravitational mass a + P is strictly negative. 

2. The field equations 

We consider a spherical shell in vacuum. The shell “separates” a portion 
of Minkowski space (“interior”) from Schwarzschild space of arbitrarily 
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fixed mass parameter M (“exterior”). The field equations relate the jump 
of the exterior curvature with the matter content of the shell, i.e. a surface 
distribution of stress-energy. It can be shown iflbll that such equations can 
be written as follows: 

(2.1) <7=- 2 |(<7 + P) 


( 2 . 2 ) 


R 2 = V(R,a) 


where R(t) is the “radius” of the shell as a function of the proper time on the 
shell, a and P are the surface (rest-frame) energy-density and the surface 
pressure (due to spherical symmetry, the two eigenvalues of the stress are 
equal, i.e. the material is necessarily a “two-dimensional perfect fluid”) and 
the “effective potential” V is given by 


V(R,a ):= 


M 2 


167T 2 P 4 <7 2 


4tt 2 PV - 1 + 


M 

~R' 


The equations (12.11) and (12.21) contain three unknowns, and become a 
closed set once a further relation has been given. This has been achieved 
in various ways in the literature; for instance, one can assign the pressure 
profile as a “a priori” known function ifTTIl and calculate the corresponding 
density profile. However, in this way, the resulting “equation of state”, i.e. 
the relationship between density and pressure which holds during the dy¬ 
namics may turn out to be unphysical. Here, we are going to consider the 
case of a closed equation of state of the barotropic type, namely P = T(cr). 
In this case, as we shall see, the equations above decouple and the dynamics 
can be qualitatively analyzed. 

The function T has to satisfy a certain set of assumptions which assure the 
physical reasonability of the model. It seems that a reasonable compromise 
between generality and physical reasonableness can be obtained assuming 
that T(cr) is a continuous, monotonic function, differentiable except perhaps 
at a = 0, where its value vanishes. These assumptions comprise a wide 
range of equations of state and, in particular, the linear-barotropic (P = T 0 er 
with a constant T 0 ) that has been widely studied in the literature lfl8ll . fTTl . 
and the polytropic case P = T 0 cr 1 ' (both these e.o.s. reduce to the newtonian 
poly tropes in the weak field regime fl9l - lt20ll ). 

If the weak energy condition is required, then solutions have to be se¬ 
lected in such a way that a > 0, and on such solutions the effective gravita¬ 
tional mass 


Ai(cr) := cr + T(cr) 


must be non negative. Therefore the requirement T(cr) > —a has to be 
added. However, shells with positive density but negative effective mass 
are of interest on their own, for instance as toy models for sources of dark 
energy (“phantom” fields, see e.g. llT2l . exotic matter and wormholes, see 
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G2), thus we shall consider also this possibility in what follows. 


3. The structure of the solutions 

We take initial data a(R 0 ) = cr 0 > 0 and consider the behavior of the 
solutions of the ordinary differential equation 

(3.1) & = 

First of all, we notice that it can be reduced to a quadrature 

r a ( R ) da 

(3.2) / —- = -21og( J R/i? 0 ). 

J* 0 F-O) 

Thus, if the function 1 /n(a) is integrable in neighborhood of +oo, the 
density diverges then this occurs at a finite non zero value R s (with R 0 > 
R s > 0) (“Case A”). If instead the function l//i(cr) is not integrable in 
neighborhood of +oo, the density can diverge only as R approaches zero 
(“Case B”). Further, a straightforward analysis shows that the following 
cases have to be distinguished: 

Case I): /r(cr) > 0 for any a: 

In this case each solution is strictly decreasing from +oo to zero (one 
among behaviors A or B occurs). Weak energy condition is always satis¬ 
fied. Since a = 0 is a solution, if dT/da is continuous in a = 0 the standard 
uniqueness theorem for ODE gives a —> 0 as R — >+oo. For completeness 
we notice, however, that one could also conceive models in which dT / da 
has a finite jump in a = 0. In this case a loss of uniqueness occurs at a 
finite value Ri (at which the solution goes to zero and can be prolonged to 
be zero, thus giving a sort of “evaporation” of the shell). We thus have four 
qualitatively different cases reported in figures 1-4. 

Case II): There exist a zero (at a p ) of the function fi(a): 

Clearly a = a p is a regular solution at which no loss of uniqueness can 
occur. Thus the behavior depends on the initial datum for a: 

(Ha) If a 0 > a p then all solutions are strictly decreasing and approach 
from above a p as R —> +oo; one among cases A and B occurs. Weak 
energy condition is satisfied everywhere. 

(lib) If a 0 < a p then all solutions are strictly increasing and approach 
a p from below as R —> +oo; their behavior “in the past” intersects a = 0 
where a loss of uniqueness occurs. Weak energy condition is violated every¬ 
where because the density is positive but /i is strictly negative. Two behav¬ 
iors must still be distinguished: if the derivative a" = —a'/R{2> + dT/da) 
has no zeroes, the solutions come out from zero with infinite tangent, if it 
has a zero the solutions have a flex point, and the matching at a = 0 is 
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smooth in this case. All in all we have the four qualitatively different cases 
reported in figures 5-8. 


4. Discussion and conclusion 

Once the behavior of the energy function a is known the corresponding 
qualitative behavior of the shell radius R(t) can be read off from the “effec¬ 
tive potential” equation (12.2k 

First of all, we notice that, due to the presence of the M/R + Ait 2 R 2 a 2 
term, the potential always diverges at the left end side of the interval of def¬ 
inition, in all the eight qualitatively different cases discussed above. The 
behavior at the other extreme is divergent as well in the cases I without 
loss of uniqueness and in all cases II (where V diverges as R goes to in¬ 
finity). Thus in all such cases the potential is a continuous curve positively 
diverging at the extremes, and it follows that it has an absolute minimum. 
The dynamics is allowed in both the neighborhoods of the extremes, and it 
is allowed everywhere if V is positive at the minimum. Bouncing or also 
oscillating behaviors can appear as well, if the minimum is negative and 
therefore there is at least one region in which the dynamics is not allowed. 

Shells with general equation of state can thus be viewed as very useful 
“toy” models for astrophysical objects, in order to investigate the effect 
of the equation of state on the dynamics in several physically interesting 
scenarios. Work in this direction is in progress. 
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Figure 3. a(R): Case IA when ^ is discontinuous in a = 0. 



Figure 4. <j(R): Case IB when ^ is discontinuous in a = 0. 
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Figure 7. a(R): Case lib when a" has no zeroes. 
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Figure 8. a(R): Case lib when o" has one zero. 
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